Detecting T-cell clonal expansions and quantifying clone survival using deep profiling of immune repertoires.

An individual's T-cell repertoire constantly changes under the influence of external and internal factors. Cells that do not receive a stimulatory signal die, while those that encounter and recognize a pathogen or receive a co-stimulatory signal divide, resulting in clonal expansions. T-cell clones can be traced by monitoring the presence of their unique T-cell receptor (TCR) sequence, which is assembled de novo through a process known as V(D)J rearrangement. Tracking T cells can provide valuable insights into the survival of cells after hematopoietic stem cell transplantation (HSCT) or cancer treatment response and can indicate the induction of protective immunity by vaccination. In this study, we report a bioinformatic method for quantifying the T-cell repertoire dynamics from TCR sequencing data. We demonstrate its utility by measuring the T-cell repertoire stability in healthy donors, by quantifying the effect of donor lymphocyte infusion (DLI), and by tracking the fate of the different T-cell subsets in HSCT patients and the expansion of pathogen-specific clones in vaccinated individuals.

Convergence, plasticity, and tissue residence of regulatory T cell response via TCR repertoire prism.

Suppressive function of regulatory T cells (Treg) is dependent on signaling of their antigen receptors triggered by cognate self, dietary, or microbial peptides presented on MHC II. However, it remains largely unknown whether distinct or shared repertoires of Treg TCRs are mobilized in response to different challenges in the same tissue or the same challenge in different tissues. Here we use a fixed TCRß chain FoxP3-GFP mouse model to analyze conventional (eCD4) and regulatory (eTreg) effector TCRα repertoires in response to six distinct antigenic challenges to the lung and skin. This model shows highly 'digital' repertoire behavior with easy-to-track challenge-specific TCRα CDR3 clusters. For both eCD4 and eTreg subsets, we observe challenge-specific clonal expansions yielding homologous TCRα clusters within and across animals and exposure sites, which are also reflected in the draining lymph nodes but not systemically. Some CDR3 clusters are shared across cancer challenges, suggesting a response to common tumor-associated antigens. For most challenges, eCD4 and eTreg clonal response does not overlap. Such overlap is exclusively observed at the sites of certain tumor challenges, and not systematically, suggesting transient and local tumor-induced eCD4=>eTreg plasticity. This transition includes a dominant tumor-responding eCD4 CDR3 motif, as well as characteristic iNKT TCRα CDR3. In addition, we examine the homeostatic tissue residency of clonal eTreg populations by excluding the site of challenge from our analysis. We demonstrate that distinct CDR3 motifs are characteristic of eTreg cells residing in particular lymphatic tissues, regardless of the challenge. This observation reveals the tissue-resident, antigen-specific clonal Treg populations.

Large-scale template-based structural modeling of T-cell receptors with known antigen specificity reveals complementarity features.

Introduction: T-cell receptor (TCR) recognition of foreign peptides presented by the major histocompatibility complex (MHC) initiates the adaptive immune response against pathogens. While a large number of TCR sequences specific to different antigenic peptides are known to date, the structural data describing the conformation and contacting residues for TCR-peptide-MHC complexes is relatively limited. In the present study we aim to extend and analyze the set of available structures by performing highly accurate template-based modeling of these complexes using TCR sequences with known specificity. Methods: Identification of CDR3 sequences and their further clustering, based on available spatial structures, V- and J-genes of corresponding T-cell receptors, and epitopes, was performed using the VDJdb database. Modeling of the selected CDR3 loops was conducted using a stepwise introduction of single amino acid substitutions to the template PDB structures, followed by optimization of the TCR-peptide-MHC contacting interface using the Rosetta package applications. Statistical analysis and recursive feature elimination procedures were carried out on computed energy values and properties of contacting amino acid residues between CDR3 loops and peptides, using R. Results: Using the set of 29 complex templates (including a template with SARS-CoV-2 antigen) and 732 specificity records, we built a database of 1585 model structures carrying substitutions in either TCRα or TCRß chains with some models representing the result of different mutation pathways for the same final structure. This database allowed us to analyze features of amino acid contacts in TCR - peptide interfaces that govern antigen recognition preferences and interpret these interactions in terms of physicochemical properties of interacting residues. Conclusion: Our results provide a methodology for creating high-quality TCR-peptide-MHC models for antigens of interest that can be utilized to predict TCR specificity.

Defective peripheral B cell selection in common variable immune deficiency patients with autoimmune manifestations.

Common variable immune deficiency (CVID) is a heterogeneous disorder characterized by recurrent infections, low levels of serum immunoglobulins, and impaired vaccine responses. Autoimmune manifestations are common, but B cell central and peripheral selection mechanisms in CVID are incompletely understood. Here, we find that receptor editing, a measure of central tolerance, is increased in transitional B cells from CVID patients and that these cells have a higher immunoglobulin κ:λ ratio in CVID patients with autoimmune manifestations than in those with infection only. Contrariwise, the selection pressure in the germinal center on CD27bright memory B cells is decreased in CVID patients with autoimmune manifestations. Finally, functionally, T cell-dependent activation showed that naive B cells in CVID patients are badly equipped for activation and induction of mismatch repair genes. We conclude that central tolerance is functional whereas peripheral selection is defective in CVID patients with autoimmune manifestations, which could underpin the development of autoimmunity.

Human thymopoiesis produces polyspecific CD8+ α/ß T cells responding to multiple viral antigens.

T-cell receptors (TCRs) are formed by stochastic gene rearrangements, theoretically generating >1019 sequences. They are selected during thymopoiesis, which releases a repertoire of about 108 unique TCRs per individual. How evolution shaped a process that produces TCRs that can effectively handle a countless and evolving set of infectious agents is a central question of immunology. The paradigm is that a diverse enough repertoire of TCRs should always provide a proper, though rare, specificity for any given need. Expansion of such rare T cells would provide enough fighters for an effective immune response and enough antigen-experienced cells for memory. We show here that human thymopoiesis releases a large population of clustered CD8+ T cells harboring α/ß paired TCRs that (i) have high generation probabilities and (ii) a preferential usage of some V and J genes, (iii) which CDR3 are shared between individuals, and (iv) can each bind and be activated by multiple unrelated viral peptides, notably from EBV, CMV, and influenza. These polyspecific T cells may represent a first line of defense that is mobilized in response to infections before a more specific response subsequently ensures viral elimination. Our results support an evolutionary selection of polyspecific α/ß TCRs for broad antiviral responses and heterologous immunity.

TCR repertoire profiling revealed antigen-driven CD8+ T cell clonal groups shared in synovial fluid of patients with spondyloarthritis.

Spondyloarthritis (SpA) comprises a number of inflammatory rheumatic diseases with overlapping clinical manifestations. Strong association with several HLA-I alleles and T cell infiltration into an inflamed joint suggest involvement of T cells in SpA pathogenesis. In this study, we performed high-throughput T cell repertoire profiling of synovial fluid (SF) and peripheral blood (PB) samples collected from a large cohort of SpA patients. We showed that synovial fluid is enriched with expanded T cell clones that are shared between patients with similar HLA genotypes and persist during recurrent synovitis. Using an algorithm for identification of TCRs involved in immune response we discovered several antigen-driven CD8+ clonal groups associated with risk HLA-B*27 or HLA-B*38 alleles. We further show that these clonal groups were enriched in SF and had higher frequency in PB of SpA patients vs healthy donors, implying their relevance to SpA pathogenesis. Several of the groups were shared among patients with different SpAs that suggests a common immunopathological mechanism of the diseases. In summary, our results provide evidence for the role of specific CD8+ T cell clones in pathogenesis of SpA.

Memory persistence and differentiation into antibody-secreting cells accompanied by positive selection in longitudinal BCR repertoires.

The stability and plasticity of B cell-mediated immune memory ensures the ability to respond to the repeated challenges. We have analyzed the longitudinal dynamics of immunoglobulin heavy chain repertoires from memory B cells, plasmablasts, and plasma cells from the peripheral blood of generally healthy volunteers. We reveal a high degree of clonal persistence in individual memory B cell subsets, with inter-individual convergence in memory and antibody-secreting cells (ASCs). ASC clonotypes demonstrate clonal relatedness to memory B cells, and are transient in peripheral blood. We identify two clusters of expanded clonal lineages with differing prevalence of memory B cells, isotypes, and persistence. Phylogenetic analysis revealed signs of reactivation of persisting memory B cell-enriched clonal lineages, accompanied by new rounds of affinity maturation during proliferation and differentiation into ASCs. Negative selection contributes to both persisting and reactivated lineages, preserving the functionality and specificity of B cell receptors (BCRs) to protect against current and future pathogens.

Subsequent malaria enhances virus-specific T cell immunity in SIV-infected Chinese rhesus macaques.

BACKGROUND: Coinfection with HIV and Plasmodium parasites is fairly common, but the sequence of infection with these two pathogens and their impact on disease progression are poorly understood. METHODS: A Chinese rhesus macaque HIV and Plasmodium coinfection model was established to compare the impact of pre-existing and subsequent malaria on the progression of SIV infection. RESULTS: We found that a pre-existing malaria caused animals to produce a greater number of CD4+CCR5+ T cells for SIV replication, resulting in higher viral loads. Conversely, subsequent malaria induced a substantially larger proportion of CD4+CD28highCD95high central memory T cells and a stronger SIV-specific T cell response, maintained the repertoire diversity of SIV-specific T cell receptors, and generated new SIV-specific T cell clonotypes to trace SIV antigenic variation, resulting in improved survival of SIV-infected animals. CONCLUSION: The complex outcomes of this study may have important implications for research on human HIV and malaria coinfection. The infection order of the two pathogens (HIV and malaria parasites) should be emphasized. Video abstract.

Accounting for B-cell Behavior and Sampling Bias Predicts Anti-PD-L1 Response in Bladder Cancer.

Cancer immunotherapy is predominantly based on T cell-centric approaches. At the same time, the adaptive immune response in the tumor environment also includes clonally produced immunoglobulins and clonal effector/memory B cells that participate in antigen-specific decisions through their interactions with T cells. Here, we investigated the role of infiltrating B cells in bladder cancer via patient dataset analysis of intratumoral immunoglobulin repertoires. We showed that the IgG1/IgA ratio is a prognostic indicator for several subtypes of bladder cancer and for the whole IMVigor210 anti-PD-L1 immunotherapy study cohort. A high IgG1/IgA ratio associated with the prominence of a cytotoxic gene signature, T-cell receptor signaling, and IL21-mediated signaling. Immunoglobulin repertoire analysis indicated that effector B-cell function, rather than clonally produced antibodies, was involved in antitumor responses. From the T-cell side, we normalized a cytotoxic signature against the extent of immune cell infiltration to neutralize the artificial sampling-based variability in immune gene expression. Resulting metrics reflected proportion of cytotoxic cells among tumor-infiltrating immune cells and improved prediction of anti-PD-L1 responses. At the same time, the IgG1/IgA ratio remained an independent prognostic factor. Integration of the B-cell, natural killer cell, and T-cell signatures allowed for the most accurate prediction of anti-PD-L1 therapy responses. On the basis of these findings, we developed a predictor called PRedIctive MolecUlar Signature (PRIMUS), which outperformed PD-L1 expression scores and known gene signatures. Overall, PRIMUS allows for reliable identification of responders among patients with muscle-invasive urothelial carcinoma, including the subcohort with the low-infiltrated "desert" tumor phenotype.

HLA variants have different preferences to present proteins with specific molecular functions which are complemented in frequent haplotypes.

Human leukocyte antigen (HLA) genes are the most polymorphic loci in the human genome and code for proteins that play a key role in guiding adaptive immune responses by presenting foreign and self peptides (ligands) to T cells. Each person carries up to 6 HLA class I variants (maternal and paternal copies of HLA-A, HLA-B and HLA-C genes) and also multiple HLA class II variants, which cumulatively define the landscape of peptides presented to T cells. Each HLA variant has its own repertoire of presented peptides with a certain sequence motif which is mainly defined by peptide anchor residues (typically the second and the last positions for HLA class I ligands) forming key interactions with the peptide-binding groove of HLA. In this study, we aimed to characterize HLA binding preferences in terms of molecular functions of presented proteins. To focus on the ligand presentation bias introduced specifically by HLA-peptide interaction we performed large-scale in silico predictions of binding of all peptides from human proteome for a wide range of HLA variants and established which functions are characteristic for proteins that are more or less preferentially presented by different HLA variants using statistical calculations and gene ontology (GO) analysis. We demonstrated marked distinctions between HLA variants in molecular functions of preferentially presented proteins (e.g. some HLA variants preferentially present membrane and receptor proteins, while others - ribosomal and DNA-binding proteins) and reduced presentation of extracellular matrix and collagen proteins by the majority of HLA variants. To explain these observations we demonstrated that HLA preferentially presents proteins enriched in amino acids which are required as anchor residues for the particular HLA variant. Our observations can be extrapolated to explain the protective effect of certain HLA alleles in infectious diseases, and we hypothesize that they can also explain susceptibility to certain autoimmune diseases and cancers. We demonstrate that these differences lead to differential presentation of HIV, influenza virus, SARS-CoV-1 and SARS-CoV-2 proteins by various HLA alleles. Taking into consideration that HLA alleles are inherited in haplotypes, we hypothesized that haplotypes composed of a combination of HLA variants with different presentation preferences should be more advantageous as they allow presenting a larger repertoire of peptides and avoiding holes in immunopeptidome. Indeed, we demonstrated that HLA-A/HLA-B and HLA-A/HLA-C haplotypes which have a high frequency in the human population are comprised of HLA variants that are more distinct in terms of functions of preferentially presented proteins than the control pairs.

New Regions With Molecular Alterations in a Rare Case of Insulinomatosis: Case Report With Literature Review.

Insulinomatosis is characterized by monohormonality of multiple macro-tumors and micro-tumors that arise synchronously and metachronously in all regions of the pancreas, and often recurring hypoglycemia. One of the main characteristics of insulinomatosis is the presence of insulin-expressing monohormonal endocrine cell clusters that are exclusively composed of proliferating insulin-positive cells, are less than 1 mm in size, and show solid islet-like structure. It is presumed that insulinomatosis affects the entire population of ß-cells. With regards to molecular genetics, this phenomenon is not related to mutation in MEN1 gene and is more similar to sporadic benign insulinomas, however, at the moment molecular genetics of this disease remains poorly investigated. NGS sequencing was performed with a panel of 409 cancer-related genes. Results of sequencing were analyzed by bioinformatic algorithms for detecting point mutations and copy number variations. DNA copy number variations were detected that harbor a large number of genes in insulinoma and fewer genes in micro-tumors. qPCR was used to confirm copy number variations at ATRX, FOXL2, IRS2 and CEBPA genes. Copy number alterations involving FOXL2, IRS2, CEBPA and ATRX genes were observed in insulinoma as well as in micro-tumors samples, suggesting that alterations of these genes may promote malignization in the ß-cells population.

HLA binding of self-peptides is biased towards proteins with specific molecular functions.

Human leukocyte antigen (HLA) is highly polymorphic and plays a key role in guiding adaptive immune responses by presenting foreign and self peptides to T cells. Each HLA variant selects a minor fraction of peptides that match a certain motif required for optimal interaction with the peptide-binding groove. These restriction rules define the landscape of peptides presented to T cells. Given these limitations, one might suggest that the choice of peptides presented by HLA is non-random and there is preferential presentation of an array of peptides that is optimal for distinguishing self and foreign proteins. In this study we explore these preferences with a comparative analysis of self peptides enriched and depleted in HLA ligands. We show that HLAs exhibit preferences towards presenting peptides from certain proteins while disfavoring others with specific functions, and highlight differences between various HLA genes and alleles in those preferences. We link those differences to HLA anchor residue propensities and amino acid composition of preferentially presented proteins. The set of proteins that peptides presented by a given HLA are most likely to be derived from can be used to distinguish between class I and class II HLAs and HLA alleles. Our observations can be extrapolated to explain the protective effect of certain HLA alleles in infectious diseases, and we hypothesize that they can also explain susceptibility to certain autoimmune diseases and cancers. We demonstrate that these differences lead to differential presentation of HIV, influenza virus, SARS-CoV-1 and SARS-CoV-2 proteins by various HLA alleles. Finally, we show that the reported self peptidome preferences of distinct HLA variants can be compensated by combinations of HLA-A/HLA-B and HLA-A/HLA-C alleles in frequent haplotypes.

A T cell repertoire timestamp is at the core of responsiveness to CTLA-4 blockade.

Biology of the response to anti-CTLA-4 involves the dynamics of specific T cell clones. Reasons for clinical success and failure of this treatment are still largely unknown. Here, we quantified the dynamics of the T cell receptor (TCR) repertoire, throughout 4 weeks involving treatment with anti-CTLA-4, in a syngeneic mouse model for colorectal cancer. These dynamics show an initial increase in clonality in tandem with a decrease in diversity, effects which gradually subside. Furthermore, response to treatment is tightly connected to the shared and public parts of the T cell repertoire. We were able to recognize time-dependent behaviors of specific TCR sequences and cell types and to show the response is dominated by specific motifs. We see that a single, specific time point might be useful to inform a physician of the true response to treatmentThe research further highlights the importance of temporal analyses of the immune response.

Benchmarking of T cell receptor repertoire profiling methods reveals large systematic biases.

Monitoring the T cell receptor (TCR) repertoire in health and disease can provide key insights into adaptive immune responses, but the accuracy of current TCR sequencing (TCRseq) methods is unclear. In this study, we systematically compared the results of nine commercial and academic TCRseq methods, including six rapid amplification of complementary DNA ends (RACE)-polymerase chain reaction (PCR) and three multiplex-PCR approaches, when applied to the same T cell sample. We found marked differences in accuracy and intra- and inter-method reproducibility for T cell receptor α (TRA) and T cell receptor ß (TRB) TCR chains. Most methods showed a lower ability to capture TRA than TRB diversity. Low RNA input generated non-representative repertoires. Results from the 5' RACE-PCR methods were consistent among themselves but differed from the RNA-based multiplex-PCR results. Using an in silico meta-repertoire generated from 108 replicates, we found that one genomic DNA-based method and two non-unique molecular identifier (UMI) RNA-based methods were more sensitive than UMI methods in detecting rare clonotypes, despite the better clonotype quantification accuracy of the latter.

T-cell tracking, safety, and effect of low-dose donor memory T-cell infusions after αß T cell-depleted hematopoietic stem cell transplantation.

The delayed recovery of adaptive immunity underlies transplant-related mortality (TRM) after αß T cell-depleted hematopoietic stem cell transplantation (HSCT). We tested the use of low-dose memory donor lymphocyte infusions (mDLIs) after engraftment of αß T cell-depleted grafts.A cohort of 131 pediatric patients (median age 9 years) were grafted with αß T cell-depleted products from either haplo (n = 79) or unrelated donors (n = 52). After engraftment, patients received mDLIs prepared by CD45RA depletion. Cell dose was escalated monthly from 25 × 103 to 100 × 103/kg (haplo) and from 100 × 103 to 300 × 103 /kg (MUD). In a subcohort of 16 patients, T-cell receptor (TCR) repertoire profiling with deep sequencing was used to track T-cell clones and to evaluate the contribution of mDLI to the immune repertoire.In total, 343 mDLIs were administered. The cumulative incidence (CI) of grades II and III de novo acute graft-versus-host disease (aGVHD) was 5% and 2%, respectively, and the CI of chronic graft-versus-host disease was 7%. Half of the patients with undetectable CMV-specific T cells before mDLI recovered CMV-specific T cells. TCR repertoire profiling confirmed that mDLI-derived T cells significantly contribute to the TCR repertoire up to 1 year after HSCT and include persistent, CMV-specific T-cell clones.

SARS-CoV-2 Epitopes Are Recognized by a Public and Diverse Repertoire of Human T Cell Receptors.

Understanding the hallmarks of the immune response to SARS-CoV-2 is critical for fighting the COVID-19 pandemic. We assessed antibody and T cell reactivity in convalescent COVID-19 patients and healthy donors sampled both prior to and during the pandemic. Healthy donors examined during the pandemic exhibited increased numbers of SARS-CoV-2-specific T cells, but no humoral response. Their probable exposure to the virus resulted in either asymptomatic infection without antibody secretion or activation of preexisting immunity. In convalescent patients, we observed a public and diverse T cell response to SARS-CoV-2 epitopes, revealing T cell receptor (TCR) motifs with germline-encoded features. Bulk CD4+ and CD8+ T cell responses to the spike protein were mediated by groups of homologous TCRs, some of them shared across multiple donors. Overall, our results demonstrate that the T cell response to SARS-CoV-2, including the identified set of TCRs, can serve as a useful biomarker for surveying antiviral immunity.

CD4+ T Cells Recognize Conserved Influenza A Epitopes through Shared Patterns of V-Gene Usage and Complementary Biochemical Features.

T cell recognition of peptides presented by human leukocyte antigens (HLAs) is mediated by the highly variable T cell receptor (TCR). Despite this built-in TCR variability, individuals can mount immune responses against viral epitopes by using identical or highly related TCRs expressed on CD8+ T cells. Characterization of these TCRs has extended our understanding of the molecular mechanisms that govern the recognition of peptide-HLA. However, few examples exist for CD4+ T cells. Here, we investigate CD4+ T cell responses to the internal proteins of the influenza A virus that correlate with protective immunity. We identify five internal epitopes that are commonly recognized by CD4+ T cells in five HLA-DR1+ subjects and show conservation across viral strains and zoonotic reservoirs. TCR repertoire analysis demonstrates several shared gene usage biases underpinned by complementary biochemical features evident in a structural comparison. These epitopes are attractive targets for vaccination and other T cell therapies.